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Abstract 

The analysis of objects living on ultrametric trees, in particular the block-diagonalization of 
4— replica matrices M a/3nS , is shown to be dramatically simplified through the introduction of prop- 
erly chosen operations on those objects. These are the Replica Fourier Transforms on ultrametric 
trees. Those transformations are defined and used in the present work. 
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Resume 



On montre que l'analyse d'objets vivant sur un arbre ultrametrique, en particulier, la diagonal- 
isation par blocs d'une matrice M al3 ' jS dependant de 4— repliques, se simplifie de fagon dramatique 
si Ton introduit les operation appropriees sur ces objects. Ce sont les Transformee de Fourier de 
Repliques sur un arbre ultrametrique. Ces transformations sont definies et utilisees dans le present 
travail. 
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Spin glasses and their typical glassy phases appear to be present in a wide spectrum of domains jl|-[|3j 
The high level of complexity inherent to their structure (field 4>c,p(x) depending upon two replicas, 
propagators G al3 ' jS [x, y) depending upon four of them) has prevented so far a systematic study of 
the glassy phase with the standard tools of field theory and in particular the renormalization group 
where objects with 6 (or 8) replicas would be needed. One recent step in the process of founding 
a field theory has been the spelling out of a Dyson like equation^, [5| relating the propagators G 
and their inverses M (mass operator plus kinetic terms), a triviality in standard field theories. To 
be more precise, the obtained relationship is not directly between G and M, but between what was 
termed the kernel (F) of G and the kernel (K) of M. What makes the algebra complicated, is that 
one cannot directly work with replicas (whose number n is to be set to 0). Instead, one is led to 
use, for the observables, a replica symmetry broken representation (RSB) usually inferred from mean 
field studies (be it with R — 1, one step RSB or R — oo as in Parisi'sj?]] ansatz). As a result, the 
inversion process of the 4— replica matrix M al3nS , becomes a highly non trivial exercise with results 
fairly involved for R = oo|J], and even much more so for a generic i?|H). 

Here we would like to show that with the appropriate choice of transformations, the previously 
obtained results, considerably simplify, and in some sense, become almost transparent. 

To that effect we introduce the notion of (i) Replica Fourier Transform (RFT) and (ii) RFT on a 
(ultrametric) tree. With those definitions, given in section 1, we show that the relationship between 
(matrix) functions and their kernels is nothing but a Replica Fourier Transformation (a double trans- 
form for the Replicon component and a single one for the longitudinal anomalous component). This 
is done in section 2 and 3. Further the same RFT, used now on the eigenvalue equations immedi- 
ately (block-)diagonalizes them (section 4). Used on the equation GM = 1, the RFT yields Dyson's 
equation relating the corresponding kernels of M and G (section 5). 

It must be mentioned at the outset that this work was prompted by the yet unpublished results of 
Parisi and Sourlas|| where a Fourier transform on p— adic numbers is introduced^, which diagonalizes 
the Replicon sector components. The RFT used here has the flexibility that also allows for a block 
diagonalization of the Longitudinal Anomalous sector. 



1 Replica Fourier Transform 

It is assumed, as in Parisi[Q], that the 2— replica object (cj) a p(S)) = q a/ 3 = q t only depends upon the 
overlap (or codistance a fi /3 = i). Likewise a m— replica object will depend upon m — 1 independent 
overlaps. The RFT|^], which is a discretized version of the algebra introduced by Mezard and 
Parisi^] defines the RFT transform A of A via 

1 See also B.Grossman, ref 
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A k = ^ Pt {A t -A t - l ) (1) 

t=k 

Here the pt's are the sizes of the Parisi boxes, po = n the replica number, Pr+i = 1. Note that 
q H P — t means that a and (3 belong to the same Parisi box of size pt, but to two distinct boxes of 
size pt+i, i-e. they belong to the (relative) multiplicity 

&t=Pt- Pt+i (2) 

Sr+i = Pr+i (3) 

Objects with indices out of the range of definition (here (0, R + 1)) are taken as null. Conversely 
one has 

j 

k=0 

and the associated (relative) multiplicity in what will be termed below "resolution" (or pseudo- 
momentum) space, 

4 = --— (5) 

Pk Pk-i 

So = — (6) 

Characteristically, the convolution 

n 

y ' Aa^B-fp — Cct/3 (7) 



7 = 1 



becomes, after RFT 



If we take 



we have 



A k B k = C k (8) 



Caf3 = Saf3 (9) 



C k = 1 (10) 
A k = 1/S k (11) 
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This convolution property conveniently allows to write for example, 



tr(A) s =nJ2Sk( K A k y , (12) 

fc=0 

to be compared with 

R+l 

^(^r = n^>(A);. (13) 

a/3 j = 

One is often using a multiplicity fj,(k), related to the (relative) multiplicity 8k introduced here by 

H(k) = nS k (14) 
Ai(0) = n— = 1 (15) 

2 RFT on a tree 

Despite their usefulness the above introduced objects are not yet tailored to fit the complexity of 
situations arising in the study of functions of more than two replicas. We need to extend the RFT 
definition to the case where the two replicas (i.e. the overlap) summed over, move on a (ultrametric) 
tree in the presence of other passive overlaps. The simplest example is the 3— replica function (see 
below section 4) 

r"-" 1 = n (i6) 

where we have a (fixed) overlap r = a n j3, and the lower index, t is the cross-overlap t 

max(a D n; /3 D /x) = t 

an/3 = r (17) 

As we move fj, i.e. t along the 3— replica tree, we now have the successive multiplicities 



Pt -Pt+i 
p t - 2pt+i 
2(pt - Pt+i) 

Indeed, in the case t = r, there are two boxes of size 
one of 0). In the case t > r there are two branches 



t < r 

t = r (18) 
t > r 

r + l excluded for \i (the box p r +i of a, and the 
of the tree available for fi (afl(i = i, PI fi = r 
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or a Pi /it = r, (3 D = t). It is thus appropriate to define effective boxes in the presence of a passive, 
direct overlap r, viz. p[ r) with 

p[ r) = pt t<r 
p\ r) =1p t r<t (19) 

and the associated RFT on the 3-replica tree 

R+l 



^ = £ 

t=k 



p[ r) (4 r) - A&) (20) 



(2D 

t=o Pt 

The RFT involved in this paper will always concern cross-overlaps, i.e. lower indices. In the case 
we are on a 4— replica tree with two passive, direct overlaps, e.g. for r < s we shall use p[ r ' (see 
sections 3.2, 4.1) with 



p\ = pt t < r 
p ( t' s) = 2p t r <t<s (22) 
p(r,s) _ ^ r < s < t 

and obvious extensions for r = s or when more passive overlaps are involved. 

3 Kernels as RFT on a tree 

Consider the 4— replica matrix M al3rfS which we choose to parametrize as follows: 

(i) on the Replicon-\ike configurations of the 4— replica tree: i.e. a PI (3 = 7 nS = r 

M ap '' iS = m;. ; ; (23) 

with the lower indices 

max(a n 7, a n 5) — u, 

u, v > r + 1 

max(/3 n 7, f3 n S) = v 

(ii) The other configurations of the 4— replica tree belong exclusively to the so called Longitudinal- 
Anomalous (LA) component 
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M a/3;7<s = Mp s = A Ml' a (24) 
where max(a D 7, a PI 5, yS PI 7, /3 D <5) = t, and where it may happen, accidentally, that r = s. 

The upper indices take values 0, 1, . . . , R (for the special problem considered here, R+l = (a Pi a) 
is excluded). Lower indices take values 0, 1, ... , R + l. These two sets of variables will also be referred 
as direct overlaps and cross-overlaps, respectively. 

We now show, that the "kernels" defined in ref.^j, [|, in terms of which it was possible to invert 
the 4— replica matrix M al3 ''~' s , are nothing but appropriate RFT's on lower indices (cross-overlaps). 

3.1 The Replicon component rM 

It was recognized quite early Jill l^l that rM^I (or the corresponding component of the inverse 
i.e. of the propagator rGQ^wss obtained via some double transform of a more elementary object, 
the "kernel" K^'J (or for the propagator). The kernel, identified with the Replicon eigenvalue 
X(r;k,l) (see below), was indeed explicitly written as|l3|] 

kinetic term + A(r; k, I) = KZ'^ = 

R+l R+l 

= E h E ft " ^ " + < L Vi) - k,l>r + l (25) 

11 — k v — L 

which we recognize now as a double RFT. Note that 

(i) we have written M instead of rM. We shall see below why this does not make any difference. 

(ii) we have not used the RFT on the tree (with passive r) since here k, I > r + 1, rendering its use 
trivial (and pedantic). 

The direct relationship is then obtained by inverting the double RFT, 

U V 

= E y k E ^ W - K Zi, - K Z+i + K Zl;l+i) (26) 

k=r + l " i = r + l 

This was, unknowingly, the type of relationship between rG^ v and that first came out in the 
unravelling of the bare propagators (for the Parisi limit, ii = 00) As already mentioned, analogous 
results have also been recently obtained, through the use of p— adic theory, by Parisi and SourlasQ. 

3.2 The LA-component A M 

With the above defined RFT on the tree (|22|) we can now write 
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M " S =£^W (27) 
fe=0 p fc 

R+l 

^ r;s = E^ r;s) [ M ' r;s - M '- s i] ( 28 ) 

Corresponding equations may be written relating the propagator G r t ' a to its kernel, obtained 
by RFT, F^' s . Equation generates the LA kernel K^' a , including the limiting value r — s. 
Conversely, knowing the kernel everywhere, we can generate the mass operator aM for all values of 
r, s. 

In the Replicon configurations of the 4— replica tree we have two lower indices and therefore 

max(ii,?j) 

E ( 29 ) 

Note that pjjT'"' 1 ^ = p^ r ' m,n ("'"^ m tne interval bounded by max(u,v). From j29|), it follows that 

A M r J v = aM^ + aM t v ' t - A M T r < r , u,v>r+l (30) 

where aMI' t is given by (27) taken at r = s. This relationship explains why on the righthand side 
of ( |25| ) one may use, indifferently, M or R M: indeed aM^I where u, v > r + 1, is shown in 
depend upon a sing/^] lower index at a time, it is thus projected out under a double RFT. 

The equations given above for the kernels (|25| , p8|) can be taken as defining the kernels once the 
primary functions M are known (e.g. by loop expansion). It is shown below that the same objects 
then block-diagonalize the eigenvalue equations (section 4) and the Dyson equation (section 5). 



4 Eigenvalue equations 

We now turn to the eigenvalue equations for the matrix M°^ ;1< . The three eigenvector classes, as 
introduced by de Almeida and ThoulessQ, write f a/3 , f al3 ^, f a0 '^ v for the L, A and R sectors 
respectively. 

4.1 The L-sector 

i M^'^f 5 = XlT 13 , (31) 

'yS 

here the sum is over all the configurations of the 4— replica tree. Spelling it out we get 

2 Note that if one wishes to attach a second lower index to G r t ' s (e.g. for the sake of faithfully representing a propagator 
by one pair of lines) it has to be min(r, s, t). 
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R+l R+l 



E \ c E E 5 * m s + 1 E ^^r* [ r = v (32) 

s = L u=r+l u = r + l t=0 J 

where the S^ s term is the contribution of the Replicon-like configurations of the 4— replica tree. 
In terms of RFT's, this writes 

R 

E {*S*cvv+i + } r = A^r (33) 

Here we have used 



r(r;s) _ „(r:s) 



Pi -Pt+i (34) 
and p' r,s ' as in ([22]). We have also used the fact that, for a sum carried over the full definition interval 
(0, R + l) or when in the R-sector over (r + 1, R + 1), one has 



R+l R+l 

^26tAt = ^2pt(At-A t - 1 )=A b (35) 

b b 



i. e. yielding the RFT at its lower bound & = 0orfe = r + l respectively. Hence in eq. (|32j) the sum 
over 8 U , S v is a double RFT. Therefore it projects out aM^ v and yields the RFT at its lower bound 
value b — r + 1 for hM£{J. For the single <5{ r ' 3 ' sum, we recover the single RFT of M t r,s at its lower 
bound 6 = 0. 



4.2 The A-sector 

I ^ M a ^ s f s ^ = Aa/ q/3;m (36) 

The idea is now the following: just like before we obtained simple relationships by taking RFT 
over lower indices, i.e. cross-overlaps, we now take one more step and do RFT on the convolution 
product of cross-overlaps. More precisely, we have the concatenations 

a/3; yS "S; " and a/3; 

or 

a/3; 7<5 yS; fi and a/3; fx 
and two other possibilities depending on the configurations of the 5— replica tree. The replica 
pairs with alike dots are those whose (cross-)overlap is to be Replica Fourier transformed. We sum 
first over the 5— replicas, at a n j3 = r, 7 n 8 = s, fixed and passive, to extract the k RFT component 
(as in going from (7) to (8) ) and then we sum sum over s to obtain 
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(37) 



Note that 

(i) the summation over the concatened lower indices yields the product K^fk (as in ^|). 

(ii) in the Replicon subspace, there is a second lower index which is freely summed over its complete 
range (r + 1, R + 1) hence yielding the b — r + 1 component of K r ' r (as in (|35|)) 

(iii) the s upper index summation, at fixed (RFT) cross-overlap k, comes out as <5s fc-1 \ 

(iv) the "monochromatic" RFT /£ corresponds to an eigevector (see (pL9|-f2 1|)) 







1 rr 



t < A; — 1 



t = k-l 



fk t>k-l 



(38) 



that is mtZZ for the cross-overlaps smaller than the resolution k, and independent of the cross- 
overlap when larger than k. 

(v) The L case identifies with k = and 



J — — 70 
po 

(vi) The full multiplicity associated with a given fk is 



1 1 



(39) 



(40) 



/x(fc) = nSk ~ n 

\Pk pk-i 

4.3 The R-sector 

i ^ M "P\i& psw = \ R f a ^» a n /3 = 7 n S 

If n v 7^ a n f3, then there is a single (independent) cross-overlap, i.e. those eigenvectors belong 
to the LA subspace. We thus have necessarily 



(41) 



an p = j n 5 = 

and we now have a double set of concatenations, e.g. 



(42) 



a/3; jS 7(5; (lis and a/3; fiu 
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where one set is exhibited and the other is its complement. Altogether four double sets depending 
on configurations of the 6— replica tree (restricted to its Replicon sector through (|4l|)). The double 
RFT associated (block)-diagonalizes the eigenvalue equations (^) (the blocks are here of dimension 
lxl, instead of (R + 1) x (R + 1) in the LA sector): 



K k-Jfkl = -Wfc 



(43) 



KX'[ = kinetic term + A(r; fc, I), k,l>r + l (44) 

The associated multiplicities can be written in term of 8 k the (relative) multiplicity in pseudo- 
momentum or resolution space, with the proviso that, at the lower bound b of the interval of definition, 
we have 



S b = - 

Pb 



(45) 



b = as for (g), for the LA sector, and b = r + 1 for the R sector of direct overlap r. Again, using 



instead of p k would amount to some trivial change in the 8 k since fe, I > r + 1. One finds 



H(r;k, I) = ^S k Si8 r (k,l) 



(46) 



where S r (k, I) = p r — (a + l)p r+1 and a = 0, 1, 2 is the occupation of boxes p r +i by the cross-overlaps 
fe, I. For further use, this is conveniently separated as 



(47) 



^ rC g(r; fc, I) = —SkSlSr 



(48) 



Hsing(r;k,l) = < 



k,l > r + 1 



fc = r + 1, Z > r + 1 



T- + 1 



The total degeneracy becomes 



k,l = r + 1 



(49) 



H K+l 



X! X! X! ^g{r;k,l) = 



71(71 — 1) 



r=0 k=r+l l = r + l 
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R R+1 R+1 

Y Msm g (r;M) = -n(R+ 1) (50) 

r=0 fc = r + l ( = r + l 

and for the LA sector, from the (R+1) x (R+ 1) blocks 

R+l 

(R+l)'^2 f j,(k)=+n(R+l) (51) 

fe=0 

yielding back the appropriate count. 

5 Block-diagonalization of the Dyson equation 

Whatever has been found above for the matrix M and the associated kernels K in the R- 

sector, K r ^ B in the LA one) can be repeated for the inverse matrix G (the propagator matrix) and 
the associated kernels F. Expresssing their relationship via 

Y^M^CP^ =5%.^ (52) 

one is now able to block-diagonalize this matrix equation. To do this, one uses 

(i) a double RFT in the Replicon sector yielding 

KZFZ = 1 (53) 
a result which can also be obtained via p— adic theory according to Parisi and Sourlasj6j. 

(ii) a single RFT in the LA sector, viz 



^(OfcM + ^r^) ^ATm+I^- 1 ') =S%. (54) 



Here 



Afc(r) 



CVV+i if k<r + l 
KZi* if k>r + l 



and in ( f54|) we have used ( |53| ) to obtain 1/Afe(r). 

After division by Si k ~^ /4 one gets the Dyson's equation, relating the kernels of inverses 

1 1 R 1 

fe Afc(r) fc A fe (s ^ Afc(r) fc 4 * 
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or 



*f = (56) 



4A*(t)' 
t=o 



if 



F7 = ~Ak(r)F^ s A k (s) (57) 

Note that if one writes out explicitly trMG, one has to sum over multiplicities, and an exact 
cancellation occurs then between the terms in 



Msing(r;M)^ ; T^M 
coming from the R-sector, and the terms 

originating from the LA-sector. 



6 Conclusion 

Assuming that via perturbation expansion we have computed M up to some loop order then, to 
obtain the corresponding G one would do the following: 

(i) compute the kernels associated with M: via (0) and Kl' s via ( pj|) 

(ii) obtain the corresponding kernel associated with G, the inverse of M, viz i*£.' ; r = 1/K^ via j53| ) 
and F^' a via a solution of (^). This last step has been shown to be analytically feasible^] if, 
as it turns out to be at the zero loop level, K r ^ s = Kk{mmr, s) 

(iii) knowing F^ one obtains rG^ v via ( p(| transposed for F,G). The knowledge of F^' a gives 
aGuyui u,v > r + 1, via ( pj| , transposed for F, G) and the other components of aG^' b via ( p7| 
transposed for F, G) . 

The sum over multiplicities in "resolution" space generated by the inverse RFT directly constructs 
/i reg (r; k, I). The cancellations described just above, that occur between R contributions with /i s in g 
weight and the (Replicon-like) contributions arising as diagonal terms in the block-diagonalized LA 
sector are already taken care of since, via the RFT on the tree, we directly (block) diagonalize the 
eigenvalue equation and the Dyson equation. 

To conclude, let us emphasize that the "conservation law"^] on pseudo-momentum (k, or k and 
I in the double RTF) which allows for the "mass-operator" diagonalization (just like the ordinary 



3 On the ultrametric tree! I.e. resulting in the pseudo-momentum equality for a 2— point function, but ultrametric 
inequalities for a p— point function, p > 3. 
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Fourier Transform under translational invariance) should play a central role in the derivation of the 
much wanted "Feynman Rules" for the field theory of the spin glass. 



The authors are grateful to G.Parisi and N.Sourlas for making available their unpublished results. 
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